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\^ ' Phenomenological implications of a minimal extension to the Standard Model are considered, 



in which a Nambu-Goldstone boson emerges from the spontaneous breaking of a global U{1) sym- 



''^ , metry. This is felt only by a scalar field which is a singlet under all Standard Model symmetries, 

c^ ■ 

f— ^ ■ and possibly by neutrinos. Mixing between the Standard Model Higgs boson field and the new 

^^ ' singlet field may lead to predominantly invisible Higgs boson decays. The "natural" region in the 

Higgs boson mass spectrum is determined, where this minimally extended Standard Model is a 

•r-t , 

^ . valid theory up to a high scale related with the smallness of neutrino masses. Surprisingly, this 

j^ ' region may coincide with low visibility of all Higgs bosons at the LHC. Monte-Carlo simulation 

studies of this "nightmare" situation are performed and strategies to search for such Higgs boson 
to invisible (Nambu-Goldstone boson) decays are discussed. It is possible to improve the signal- 
to-background ratio by looking at the distribution of either the total transverse momentum of 
the leptons and the ^t, or by looking at the distribution of the azimuthal angle between the "^t 
and the momentum of the lepton pair for the Z- and Higgs-boson associated production. We also 
study variations of the model with non-Abelian symmetries and present approximate formulae for 
Higgs boson decay rates. Searching for Higgs bosons in such a scenario at the LHC would most 
likely be solely based on Higgs to "invisible" decays. 



1 Introduction 

In a seminal paper published in 1962, Goldstone, Salam and Weinberg [1] proved that the 
physical particle spectrum of a theory in which a continuous, global symmetry is sponta- 
neously broken must contain one massless, spin-zero particle for each broken symmetry. 
Massless particles of this type, today called Nambu- Goldstone bosons (NGB), were first 
theoretically discovered in particular models by Goldstone [2] and Nambu [3] . In the follow- 
ing, they will collectively be denoted by the symbol J . NGBs have the peculiar property 
that they couple to the divergence of the current 2^{x) associated with the symmetry that 
is broken. This coupling has a strength which is inversely proportional to the scale of 
symmetry breaking F, 

Ant = ^J{x)-d,3'{x). (1.1) 

This form of interaction is invariant under the shift transformation, J ^ J + uj, where uo is 
an angle that parameterizes different vacuum field configurations. Since NGBs typically are 
amongst the lightest particles in a theory a large fraction of the other particles can decay 
into them through eq. fll.ll) . For this decay to occur, these other particles, possibly scalars 
and/or quarks and leptons, must be charged under the same spontaneously broken global 
symmetry. It was first proposed by Suzuki and Schrock [4] that if the Standard Model 
(SM) Higgs boson mixes with such a new scalar particle then it must have a decay channel 
into a pair of NGBs {J J) provided that the scale F is of the order of the electroweak 
gauge boson masses, F ^ 100 GeV. If such a Higgs boson decay exists then it should be 
searched for at colliders. 

The basic idea underlying this article is the existence of an additional global "phantom" 
symmetry, Gp = U{l)p (P stands for "Phantom"), that is spontaneously broken at some 
scale F. Then following eq. fll.l|) . J' will couple to all fermions (/) that are charged under 
Gp since d^j'^ = mff'j^f. This coupling will be proportional to nif/F. In the literature, 
there are three famous types of Nambu- Goldstone bosons: axions [5], familons [6] and 
majorons [7] and their associated broken symmetries are the Peccei-Quinn symmetry [8], 
and the family and lepton number symmetry, respectively. In the former two cases the 
global symmetry is carried by both quarks and leptons and in the latter case by leptons 
only. However, considerations of energy loss in stars, supernovae and/or in terrestrial 
collider experiments [9] conclude that F > 10^^^'^ GeV in these popular cases. This bound 
constrains the decays of Higgs particles into the NGBs of the aforementioned models to 
be completely unobservable at colliders. Recently, a Majoron model has been considered 
where lepton number is spontaneously broken at the electroweak scale but in accordance 



with astrophysical bounds [10], however we will not consider this class of models here. 
A different situation arises if we assume that such additional NGBs, if existent, must 
exclusively couple to phantom (SM gauge singlet sector) fields. 

It is important to note that the requirement of renormalizability poses some constraints 
on such a hypothetical phantom sector. In particular, it demands that the only places 
where a phantom sector can make connections to the SM are the Yukawa interactions 
of neutrinos and Higgs bosons and the H^H "mass" term. Therefore, the only relevant 
phantom sector fields are a right-handed fermion (possibly coming in three copies) and 
(in general complex) scalar fields. This immediately triggers some thought on implications 
for neutrino masses. For them, there are two possibilities: Majorana or Dirac masses. 
The Majorana see-saw mechanism [11] in fact is nothing but a type of phantom sector. 
However, as already discussed, in the simplest models the possible spontaneously broken 
global symmetry is lepton number - clearly not a purely phantom sector symmetry. So, 
what about the Dirac case? Sticking to the same principle that leads to suppressed neutrino 
masses in the Majorana see-saw scenario, an analogous non-renormalizable operator can 
be constructed. It reads 

C = (tlm^ . (1.2) 

In the model proposed in this article, some (purely phantom sector) symmetry Gp, prevents 
the interaction L ■ Hur from providing neutrinos with electroweak-scale masses. Then, 
eq. 01.21) results in acceptably small Dirac neutrino masses after spontaneous symmetry 
breaking of Gp (and the extended SM gauge group Gsm) at ($) ~ (if) ~ 100 GeV provided 
that A ~ 10^^ GeV. Here, the field H (where H = ia2H*) is the standard model SU{2)l 
Higgs doublet and "■" denotes the inner product within Gsm or Gp. A renormalizable model 
resulting in the effective operator of eq. (II. 2p was first built by Roncadelli and Wyler [12]. 
It has been recently shown in ref. [13] that this model would lead to successful baryogenesis 
via Dirac leptogenesis [13-15] if 0.1 GeV < ($) < 2 TeV. 

It is worth noting that this particular NGB evades many bounds applying to other 
species of NGB since the only fermions transforming under Gp are the z/^j, and the coupling 
between the NGB and the neutrinos is proportional to ^ ~ ■^. This is too small to affect 
neutrino flavour oscillations through u —* u + J' [16]. 

It is not unreasonable to suppose that the effects of the phantom sector may already 



have been seen Q| in experiments revealing that neutrinos have small masses. 

The existence of such a phantom sector may also be responsible for electroweak symme- 
try breaking. This has recently been emphasized by Patt and Wilczek [17] and also by the 
authors of ref. [13]. Their argument is based on the fact that no symmetry principle can 
forbid the mixing of the Higgs sector with the phantom sector through the renormalizable 
link operator 

Aink = r/iftij$t$. (1.3) 

Eq. (11. 3p suggests that the phantom sector field $ triggers spontaneous electroweak sym- 
metry breaking, i.e. {H) = v ^ 246 GeV once it develops a vacuum expectation value 
(vev), ($) = a. This holds true even in the absence of any tree-level Higgs mass term, 
fi'^H'^H [18-20]. Furthermore, it is exactly the mixing term of eq. (11. Sp that causes the Higgs 
boson to decay into a pair of NGBs, H —>■ J J . Since the J^s, interact only very weakly 
with matter, this decay effectively constitutes an invisible decay of the Higgs boson. 

Of course, this discussion could be generalized to non-Abelian groups. However, for 
simplicity here and onwards the simplest group Gp = f/(l)p is assumed. The Noether 
current associated with this symmetry is j^i^x) = ■j<l'*9^$. The phantom field $ can be 
expanded about its vev a in the usual fashion, 

$(x) = e^^(^)/"[a + 0(x)]/V2. (1.4) 

Using eq. (Il.ip the interaction between the massive Higgs boson 0(x) and the NGB is 
found to be £int = ■^<t>{d^J)'^. The scalar potential is composed of the usual quadratic and 
quartic terms for H and $ as well as the link term of eq. (II. 3p . It is independent of J i.e. 
V{H, $) = l^(/i,0), where h is the neutral field component of the SU {2) L-Kiggs doublet. 
The fields h = OnHi and cf) = 0i2Hi are rotated to their physical mass eigenstates, if,, with 
an orthogonal rotation matrix O (CP-conservation is assumed). After setting particles on 
their mass shell, £int becomes [4], 

2 
777 

Ant = -^ 0,2 H,{X) J{:X) J{X) , (1.5) 



where i = 1,2 in this minimal Gp = t/(l)p scenariqj. In this case the rank-2 matrix 



E 



^ One should also notice that, like i/s, the three right-handed neutrinos being SM-gauge singlets are the 
only light fermions that obey the shift invariance, vji -^ vji-^ uj where w is a Grasmann-type parameter. 
It may be tempting to interpret the vr^ as Goldstinos of an N f (with n/ being the number of vr flavours) 
supersymmetric phantom sector where the J7s belong to the same supermultiplet. 

^ The link term of eq. (|1.3p also gives rise to quartic HiHjJJ{i,j — 1,2) couplings which are given 

in Fig. [TT] of [Appendix A[ These couplings contribute to the decay, H2 -^ HiJJ. However, the decay 
rate for this channel is on the order of 10~^ GeV or less for benchmark scenarios considered in this paper. 
Hence, they will be completely neglected in the analysis presented here. 



O contains one mixing angle 9. Following the notation of ref. [13], it will be fixed by 
Oi2 = — O21 = sin^ and On = O22 = cos 9. The limit rj = implies 9 = 0, i.e. no 
mixing between the SM-Higgs and the phantom sector scalar fields. Obviously, in this limit 
the Standard Model is recovered. This article assumes a convention where rriHi < rnH2- 
Trading the vev of ($) = a with the more familiar tan/? = v/a, the free parameters of the 
model read 

^TT-Hi , fnH2 5 tan6' , tan/?. (1.6) 

Eq. (11.51) is the equation underlying all phenomenological analyses in this paper. It 
describes Higgs boson decays to the almost sterile NGB particles. There exists an extensive 
body of literature, which addresses various techniques for discovering an invisible Higgs 
boson at colliders. They can be assembled in three main strategies: 

• Studying the recoil of the Z-gauge boson in the associated Z + Hi production pro- 
cess. Experimental results from LEP are summarized in [21] and simulations have 
been performed in [22] . A study for this process at the Tevatron has been performed 
in refs. [23,24] with the result that the collider needs substantially more integrated lu- 
minosity to improve the current LEP exclusion limit. Parton level simulation studies 
for the LHC exist in Refs. [24-26]. Further hadron level/detector simulation studies 
for the LHC are currently under way [27,28]. 

• Vector boson fusion (VBF) processes. As suggested by Eboli and Zeppenfeld [29], 
this has now been simulated at hadron/detector level for the LHC [27,30]. 

• Central exclusive diffractive production has been studied for a particular model in 
ref. [31]. 

It should be noted that in all the above analyses only models with only one Higgs boson 
decaying completely invisibly were considered. 

In this article the focus will be put on the first two search channels, namely ZH produc- 
tion and VBF. In both cases, the coupling of the Higgs to the gauge bosons is crucial. In 
the model considered here, only the SM-like scalar field h, belonging to the SU{2)l Higgs- 
doublet, couples to vector bosons V . The corresponding SM coupling constant [(^iify y]sM 
is rescaled with the mixing angle such that 

gUiVV = Oil [gHvv]sM ■ (1-7) 

Since the matrix O is real and orthogonal, its elements are smaller than unity. This 
immediately implies that all Higgs production cross sections and/or decay rates (to SM 



particles) in this model are suppressed relative to the SM by a factor Of^. However, because 
of the orthogonality condition, X]JC)«i]^ = X]i[^«2]^ = 1- Therefore, if for example Hi is 
invisible then the other Higgs H2 tends to be visible and vice-versa. Is this a no-lose 
theorem for a Higgs boson discovery in this type of model? This is one of the questions to 
be addressed in this paper. 

More specifically, in this article the following two questions will be discussed: 

Ql: Is there any window in the parameter space (11. 6p where LEP failed to exclude both 
Higgs-bosons for niH. < 114 GeV ? 

Q2: Is there any (natural) window in the parameter space (II. 6p where both Higgs bosons 
would hide undetected at LHC? 

In this context "natural" means that the theory has a positive definite potential, with 
perturbative (non-trivial) couplings up to a high cut-off scale A ^ 10^^ GeV, where the 
mechanism for naturally light neutrino masses may be expected [recall eq. (11.21) ]. Therefore, 
we begin our analysis with Section [2] where stability and triviality bounds are analyzed and 
plotted together with electroweak p-parameter constraints. In Section[3], we answer question 
Ql. We derive analytical formulae for the Higgs boson to "visible" (7?.^) and "invisible" 
(T^) decay rates and plot predictions of the model against experimental LEP exclusion data 
for Higgs masses less than, approximately, 114 GeV. A possible scenario explaining the LEP 
Higgs boson excess is also discussed in this section. In Section HJ we extend the region of 
validity of (7?.^) and (T^) to heavier Higgs boson masses, and justify five benchmark points. 
Next, in subsections 14. 1114. 2[ we perform a detailed Monte-Carlo simulation for signals at 
these points and their backgrounds, and we discuss possible strategies useful for further 
theoretical and experimental consideration. Furthermore, in Section El extensions of the 
Abelian to non-Abelian phantom sectors and some consequences relevant for Higgs boson 
phenomenology at the LHC are discussed. A discussion of our findings together with some 



remarks for alternative scenarios is presented in Section [71 In Appendix A , we display the 
relevant Feynman rules of the Abelian model. 



2 Stability and Triviality Bounds 

In the minimal phantom model, the set of physical parameters in eq. (11.61) can be writtero 
in terms of the renormalization group running parameters Ah H'^-, rj H^ $^, and A$ "l*^: 



1 ml^ cos2 9 + ml^ sin^ 6 
1 {mjj^ — mjjj sm{26) taxi f3 



^H = -,^ ^^ , (2.1) 



2 V 



2 



:2.2) 



1 mjj_^ sin^ 6 + mj^^ cos^ 6 



A$ = -^^ :-^ tan^/3 (2.3) 

2 f^ 

with V ~ 246 GeV. Notice that in the limit where both tan f3, tan 6^0 the phantom sector 
completely decouples from the SM scalar sector. Also, note that A$ depends quadratically 
on tan/3 and the Higgs boson masses. This implies that in the case of non-zero Higgs mixing 
there is always an upper bound on tan/? if the theory is required to remain perturbative. 
For example, if tan 6' = 1 and itlh ^ 200 GeV then tan/5 < 2. In all of our plots only the 
case tan/3 = 1 is considered although, as already explained, higher values of tan/3 would 
further reduce the number of visible Higgs events. 



a. 



There are twqj classic, "theoretical" constraints on models that have been worked out 
numerous times in great detail for the SM and in many of its extensions [32]. Firstly, the 
triviality constraint is essentially the requirement that the couplings in eq. (12.11) - eq. (12. 3p 
stay perturbative up to a certain scale A^ ^ v. Secondly, the vacuum stability constraint 
demands that the potential is bound from below up to a scale Ay ^ v. Applying both 
constraints yields A7-,Av' < 10^^ GeV, where we recall the discussion following eq. (11.21) . 
The vacuum stability bound can be reduced to the requirement 

4A^(Q)A$(Q)>r7(Q)2, (2.4) 

at all scales Q < Ay- 

The running parameters are defined at the scale Qq = Mz and then evolved up to 



■^We adopt the notation of ref. [13]. 

*The unitarity constraint here is avoided by assuming that aU quartic couphngs are in a perturbative 
region, A < 1. 



higher scales with the following 1-loop renormalization group equations [33, 34] 



1 a 2 d^H 
''" dt - 


= V' 


+ 24A^ + 12AF/-6r/ 
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dt 
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12A// + 8A$-4r/ + 6ri- 


-li^9l + 9")\ 
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2 '^'^* 

167r — — = 


= 2 


7' + 20A|. 







(2.5) 

Here, t = IhQ/Qq, g' and g2 are the f/(l)y and SU{2)l gauge couplings, respectively, 
and Yt is the top quark Yukawa coupling. We ignore all other Yukawa couplings because 
their effect in the running is negligible. The equations for Y^, g' and g2 are well known [35] 
and are left out for brevity. It is worth noticing that the parameter rj is multiplicatively 
renormalized at one loop. Although there is no particular reason for r^ = 0, if this is the 
case at one energy scale then this will remain true at all energy scales. 

Fig. [T] shows the light Higgs boson mass miji vs. mH2~^T^Hi plane for tan/? = tan^ = 1 
where the background colours show the scale of new physics A required either by positivity 
or triviality (whichever is lower). The curved contour shows the 95% CL. upper limit on the 
combined Higgs boson masses from precision electroweak data (see corresponding formula 
in ref. [13]). Fig. [T] should be compared with Fig. Oof Section HI to see the correspondence 
between easily accessible regions at the LHC and regions with a potentially high effective 
theory cut-off. The light (light green) shaded parameter region of Fig. [1] is what we will 
coin the natural region throughout this paper. 

3 LEP searches 

The LEP experiments searched first for visible Higgs boson events in the Higgsstrahlung 
process e^e~ -^ ZH with the Higgs boson decaying to 6-quarks and leptons £ in final 
state stemming from the Z boson decay. They presented [36] 95% CL. upper limits for a 
parameter 71? (in their notation S95), defined as the ratio of the number of Higgs boson 
events expected in any given model to the number expected in the Standard Model for a 
Higgs boson with an identical mass, as a function of the Higgs boson mass. An important 
point to note in this context is that TZ^ only counts "visible" events. In particular, the data 
on decays to 6-quarks will be used in the following. Then the TZ^ parameter translates into 

' ~ a(e+e- ^hX) Br(/i ^ YY) ' ^ ' 

where i = 1,2, X are the remnants associated with the production of a Hi or h (the SM 
Higgs boson) and YY could be in principle either bb, or tt, but not J J . However, in the 
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Figure 1: The light Higgs boson mass tuh-^^ vs. r/i/Zj —f^Hi plane for tan f3 = tan 6 = 
1, showing the expected cut-off A of the effective theory taking the triviality and 
positivity of the potential into account (the lower of either A ^ or Ay is shown) . The 
curved line shows the 95% CL. upper limit on the Higgs boson masses stemming 
from precision electroweak data. 



framework of the particular model studied here, another possibility is that YY = HjHj. 
Exclusion limits in this case have been presented in [36]. 

The four LEP experiments [37] also performed searches for acoplanar jets (as signal for 
Z{—>- qq) H{—^ invisible) or leptons (as signal for Z{—>- ii) H{—^ invisible), with i = e,n, 
apart from the DELPHI-coUaboration which also used r's in the final state. In all cases, the 
emergence of invisible decay products of the Higgs boson is identified with the production 
of missing energy (f!). Their study resulted in an upper limit on the branching ratio of 
H —^ invisible as a function of the Higgs mass, multiplied by the production cross-section 
normalized to the rate expected from a SM Higgs decaying completely invisibly. In our 
case, this limit places constraints on the parameter 

cr(e+e- -^ HiX) 






Br H, 



JJ) 



(3.2) 



a{e+e-^hX) 

where again i = 1,2, h is the SM Higgs boson and X are the remnants associated with the 
production of Hi or h at LEP. 
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A further important constraint comes from the OPAL collaboration who performed a 
model-independent analysis of the Higgs sector at LEP [38] . They searched for the generic 
process e~^e~ -^ ZS^ where S^ is a completely neutral (and hence invisible) scalar boson. 
Since this analysis is independent of the eventual fate of the Higgs candidate it bounds the 
parameter 

(3.3) 



a[e ' e 



a{e+e~ -^ Zh) ' 
as a function of the Higgs boson mass. In this model Si = cos^ 9 and S2 = sin^ 9. 

Particularly simple expressions may be derived for TZf and T^ in the minimal phantom 
scenario provided that the narrow width approximation may be assumed and that the 
Higgs boson to off-shell gauge boson decay modes may be neglected. Our analytical findings 
closely follow the model-independent analysis of ref . [39] . Consider the case where YY = bb 
in eq. (13.11) . For simplicity let us assume that the decay H2 —>■ HiHi is kinematically 
forbidden, i.e. itihi > mH2l'^- In this case Br(ifj -^ bb) + Br(i7.j -^ J J) ~ 1. Applying 
this to eq. (13.11) in the LEP search region, mH^/2 < tjihi ^115 GeV and after some algebra 
we arrive at 



nj 



ni 



(1 777 \ 

1 + ^ tan^^tan^/?) 
12 mi / 



-1 



fl + cot' 



1 m 



'H2 



12 ml 



cot^ 9 tan^ /3 



Firstly, notice that the number of Higgs boson events where the Higgs boson decays to bb 
(or indeed any other visible mode) are always suppressed relative to the SM prediction in 
which Br(/i -^ bb) ^ 1 for this particular Higgs boson mass region. Secondly, the number 
of visible Higgs boson events decreases in the limit m^i ^ iTib. Note also that if tan/5 > 1, 
TZf receives an additional suppression. 

The importance of the Higgs boson to invisible decay and of model-independent Higgs 
boson analyses are highlighted when we consider the example tan 6' = 1 and tan/? = 2 
where we obtain 7^^ = 0.012 for m/^. = 50 GeV. In principle, this is within the region 
allowed by LEP "visible" Higgs search data [36] which excludes 0.015 < TZf < 0.2 for Higgs 
masses in the range 12 GeV < m,H, ^100 GeV. 

This could have been a "nightmare" scenario; LEP would have completely missed the 
Higgs sector! Fortunately, this nightmare is averted by both the LEP Higgs boson to invisi- 
ble searches and the OPAL model-independent Higgs boson search, because the former, for 
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Figure 2: LEP excluded regions (at 95 % C.L.) (in dark blue) for tan 6* versus 
Higgs boson masses tuhi (left plot) and 7/1^2 (right plot) for the minimal phantom 
scenario with tan/? = 1. Both searches are clearly complementary to each other 
in this scenario. 

instance, sets bounds on T^. In the relevant LEP mass region, niH^/^ < mui ^115 GeV, 






cos^e 
sin^^ 






(3.4) 
(3.5) 



Setting i?f — > implies that T^ +T2 ~ 1. LEP searches for invisible Higgs bosons exclude 
7^^ > 0.5 for masses below 110 GeV, rriH^ ^ 110 GeV. Therefore, it is unlikely that there 
are two invisible Higgs bosons in the LEP search region with masses mn^ ^110 GeV. This 
answers question Ql posed in the introduction. 

In addition, using the model-independent analysis of OPAL [38], niHi ^ 85 GeV is 
excluded for sf > 0.5. Since either sf = cos^ 9 > 0.5 or S2 = sin^ 9 > 0.5 for any given 9, 
OPAL excludes the case where both mu^ < 85 GeV and mH2 ^ 85 GeV, independently of 
how the Higgs bosons actually decay. 

It is interesting to note that one Higgs boson could still be hidden in the LEP search 
region even with these strong constraints, while the other Higgs boson then would wait for 
its discovery in the allowed region out of reach of LEP. 
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The results of a detailed analysis of this model, including visible, invisible and model- 
independent LEP bounds [36-38] are summarised in Fig. O This numerical analysis con- 
firms the analytical findings above. With tan/? = 1, a light Higgs boson (Hi) with a mass 
as low as 65 GeV could have escaped unobserved at LEP if tan 6' > 2. For the same range 
of tan 6' the other Higgs {H2) is constrained to be heavier than 114 GeV. From inspection 
of Fig. [2] we can define a LEP-allowed benchmark scenario Bl for the phantom model 
presented here, namely: 

Bl : ruH, = 68 GeV , niH^ = 114 GeV , 

tan^ = 2 , tan/? = l. (3.6) 

In this case one Higgs boson is buried, undiscovered in the LEP search region due to the 
small values of T^f and T^ which have to satisfy TZf + T^^ = cos^ 6 = 0.2 following eq. (13.41) . 
The other, heavier Higgs has TZ^ = 0.06 and T2 = 0.74. With this set of parameters, very 
few H2 events are SM-like decays into "visible" final states and instead H2 decays mainly 
into "invisible" NGBs. This scenario could well be classed as a (LEP) nightmare! 

As yet, no combined LEP limits exist on invisibly decaying Higgs bosons with masses 
below rriH = 90 GeV. Therefore, for tuh < 90 GeV the limits presented here are based on 
the individually published results from each experiment. However, some of the individual 
studies do not cover the whole Higgs mass range considered here and so the best available 
limit is used for any given Higgs mass. This is one of the causes of the sharp edges in 
Fig. [21 Clearly, a future combined LEP analysis may well exclude the benchmark Bl 
which lies close to being ruled out by ALEPH [37] which considered Higgs masses down to 
rriH = 70 GeV for which T^ ~ 0.1 is excluded. 

3.1 A digression: 2.3 a LEP Higgs search excess 

The LEP experiments established a small, 2.3cr effect in their Higgs boson searches corre- 
sponding to a Higgs boson mass of about 98 GeV [40]. Explaining this excess would require 
a value of TZf ^ 0.2, ruling out a Standard Model Higgs boson as plausible explanation. 
It is possible to provide a candidate Higgs boson in the phantom model discussed in this 
publication, which would have produced such an effect in the LEP data. Fig. [3] shows the 
allowed region in the tan^ vs. tan/5 plane for ttih^ = 98 GeV. The allowed region is tightly 
constrained because of the searches for invisible Higgs bosons at LEP in this mass region. 
At the relatively small values of tan /? still allowed, the main reason for such a small value 
of TZi is Higgs mixing rather than the extra invisible decay mode suppressing the Higgs 
branching ratio. 
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2.3 o- LEP Higgs Search Excess at mH=98 GeV 
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Figure 3: The tan^ vs. tan/3 plane for a lightest Higgs mass of 98 GeV. The 
lightest region indicates where the 2.3(7 effect in the LEP Higgs searches could be 
explained whilst still being consistent with other LEP Higgs search data (such as 
the search for invisible Higgs bosons). 

Fig. m shows the constraints on this region of parameter space coming from considering 
the triviality and positivity of the potential. For a suitably heavy mH2 ^ 210 GeV, most 
of the region suggested by the LEP excess is described by a theory which could be valid 
to scales as high as 10^^ GeV. When tan^ ~ 1, as tan/3 -^ (a — > cxd) the values of 
7] and A$ tend to 0. Looking at eqs. fl2.5p it is apparent that small values of t] and A$ 
will be relatively stable under renormalization group evolution since, for example, rj is 
multiplicatively renormalized. Higgs masses around the electroweak scale are maintained 



in this limit because n^ -^ oo whilst /x| 



However, because 77 — ;> the $ and H 



sectors are almost decoupled so that potentially destabilizing diagrams with a heavy H 
will be proportional to rj and not greatly affect the mass of ^. Of course, the model is still 
quadratically sensitive to a high cut-off scale and thus still possesses the hierarchy problem 
of the SM. 

Note that the second Higgs boson mass is restricted by the upper limit on Higgs boson 
masses from precision electroweak data [13], however for mH2 ^ 210 GeV the whole region 
suggested by the LEP excess is free from this constraint. Clearly further data would be 
required before this effect could be taken more seriously. 
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mHi = 98 GeV, m^2 = 210 GeV 
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Figure 4: The tan 6 vs. tan f3 plane for Higgs boson masses of rriHi = 98 GeV and 
mH2 = 210 GeV. The region enclosed by the black line indicates where the 2.3a 
effect in the LEP Higgs searches could be explained whilst still being consistent 
with other LEP data for Higgs boson searches. The background colours indicate 
the scale of the expected cut-off A, of the effective theory taking the triviality and 
positivity of the potential into account. Darkly shaded regions are excluded by 
LEP Higgs search data. Gontours for 7^| are shown in white. 

In the next chapter we will address the question of whether the LHG has the sensitivity 
required to discover these scenarios, in particular the potential nightmare Bl. The possible 
existence of other challenging scenarios with heavier Higgs bosons will also be examined. 

4 LHC: expectations and strategic searches 

In the LHG search region, the parameters TZ^ and T-^ can be defined by expressions similar 
to those in eqs. (13.11) and (13.21) . respectively, with the obvious replacement of the elec- 
tron/positron initial state to a proton/proton initial state and YY = 'yy,bb,VV,gg, etc.. 
Two categories for the ratios TZj may be distinguished: (a) a region where m/^. < 2 my 
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and Hi decays dominantly into hh and ( h) a region where tti^. ~ 2 my and the Hi decays 
dominantly into a gauge boson pair, VV , with V = Z^ W . 

In case (a), and under the assumption that gauge bosons are produced on-shell, analyt- 
ical approximations for TZ^ are identical to those studied in the previous chapter. On the 
other hand, assuming a common gauge boson mass my, in region (6), we obtain 

1 -1 
1 + tan^ e) (l + — -— tan^ 6 tan^ (3) 



nl 



nl 



;i + cot^ 6) (l + -^ cot^ e tan^ (3 + IM_ """^'L,. (^^^ 9 - tan/3)' 
V 3^(x2) g{x2) (l + cot2 0)2' 



(4.7) 



where Xi = ■my/m'\j^^ and g{x) = (1 — 4x + 12x^) (1 — AxY'"^. The last term in eq. (14. 7p 
is the contribution from the heavy Higgs boson decay H2 -^ HiHi [34,41]. Furthermore, 
y = m'\j^/m'\j2 and f{y) = (1 + Ay + Ay'^) (1 — AyY^"^ 0(1 — Ay). Imposing some constraints 
to this analysis (see section [2]), the mode H2 — > HiHi will not be important in further 
discussions. 

It is apparent from eq. ( 14. 7p that a certain suppression of the observable rates (Jlf) 
is always present. Its origin is twofold. Firstly, the couplings between the Hi and SM 
fields are always suppressed because of the mixing matrix O. Secondly, the decay widths 
of the Higgs bosons are enhanced by the extra decay mode Hi -^ J J . The contribution 
of this additional decay mode is increased at large tan /3 and for tan /? = 10 and tan 6* = 1 
the suppression of visible events is always more than 90% for vnu^ ^ 200 GeV. However, 
as we have already remarked in section [21 high values of tan /3 result in non-perturbative 
couplings and will therefore not be considered in this article. 

What then would be a nightmare scenario for the LHC? At present both the ATLAS 
and CMS collaborations have performed studies, at detector simulation level, to explore 
the discovery potential of their apparatus for both SM-like Higgs bosons which decay to 
visible final states, see e.g. [42], and Higgs bosons decaying to invisible final states, for 
example [28,30]. These studies are sensitive to the ratios V}^ and T^ as functions of 
the Higgs boson mass. For example, looking at the simulation results for the LHC with 
£ = 10(30) fb^ integrated luminosity we estimate (with naive scaling) that it would be 
difficult to discover a visibly-decaying Higgs if signal event rates were 30% (20%) of that 
expected in the SM (7?.^ < 0.3(0.2)). Furthermore, studies of the sensitivity of the ATLAS 
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detectoio to invisibly decaying Higgs bosons suggest that after C = 10(30) fb ^ integrated 
luminosity ATLAS could exclude Higgs bosons with T^ > 0.30(0.25) at 95% C.L. [28,30]. 

To further illustrate the necessity of the Higgs boson to invisible searches in this mini- 
mal phantom scenario, in Fig. [5] areas on the tti^i vs. mj^j — rriHi plane are plotted where 
TZf > 0.3 and/or T^ > 0.3. These limits define naive regions, where Higgs bosons will ex- 
perimentally be accessible at the LHC, either in visible or invisible search channels. Clearly 
at this stage in this study these limits are assumptions, and in fact the true experimental 
reach of the LHC will not be known until after it has been running for some time and 
predictions for the levels of backgrounds have been confirmed (or not). These assumptions 
do, however, serve as a good first estimate on which to justify the further study undertaken 
here. 

In producing Fig. [5] all Higgs boson decay modes including decays to off-shell vector 
bosons have been considered. Different colours indicate regions where either one, both 
or no Higgs bosons can be seen in different channels. It is clear that a truly challenging 
region for LHC region remains where TZj < 0.3 and T^ < 0.3. This motivates the further 
more detailed Monte Carlo analysis in the later sections of this article, which probe more 
carefully the possibility of discovering a Higgs boson in the invisible search channel when 
T^ < 0.3. 

Using equations (13. 4p and (13. 5p it is easy to see that 7^^ + 7^2 + ^^ + '^^ = 1- The 
following no-lose theorem then exists: If experiments can discover a Higgs boson over the 
whole range of 7^^ down to TZJ = 0.25 or over the whole range of T^ down to T^ = 0.25 
then at least one Higgs boson should be found. 

Without real data, estimates of the capabilities of experiments like ATLAS and CMS 
may easily be too optimistic or too pessimistic. Therefore in this publication, a constructive 
approach is taken. The phantom model has been added to the Monte Carlo event generator 
Sherpa [44], ready to be used when real data arrive. For now Fig. [5] may be used to define 
additional benchmark scenarios, some in potentially nightmarish regions, and these points 
can be studied in more detail. The particular scenarios are displayed in Table [H 

The LO branching ratios for both Higgs bosons are presented in Table[2j These ratios are 
in agreement with the analytical LO expressions in eq. (14.71) and Fig. [5l and the discussion 
following them. The most optimistic benchmark point is Bl and the most challenging one 
isB5. 

Prospects for discovering the Higgs bosons in the various benchmark scenarios B1-B5 at 
the LHC will be studied in the following. Theoretical vacuum stability and triviality bounds 



^Similar studies exist for the CMS detector [43]. 
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Ri^ > 0.3 and Tj^ > 0.3 are 'accessible' 



70 



60 



50 



> 

CD 

O 40 



30 



20 



10 








100 120 140 



1 60 1 80 200 
ahhi (GeV) 



220 240 



Figure 5: Regions in the rriHi vs. tuh-^ — f^Hi plane for tan/? = 1 and tan 6' = 1, 
where different Higgs bosons are "accessible" at LHC We define that a given H^ 
is accessible if either TZ"^ > 0.3 or 7^^ > 0.3. In the dark (blue) regions both Higgs 
bosons are accessible. In the white (beige) region no Higgs bosons are accessible. 



tan 6 = 1 tan j3 


= 1 


Benchmark 


mnAGeY) 


m^,(GeV) 


B2 


112 


130 


B3 


140 


165 


B4 


160 


190 


B5 


185 


190 



Table 1: Four LHC benchmark scenarios for the phantom model. 

as well as bounds from fitting electroweak (EW) observables have already been presented 
in Section [2l All benchmark scenarios selected in Table [1] satisfy the EW constraints and 
in some the effective theory may be valid even to scales as high as the Planck scale. 
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Benchmark 


Higgs 


TtotiGeV) 


hi 


w+w- 


ZZ 


JJ 


Bl 


Hi 
H2 


0.041 
0.051 


0.694 
3.567 


0.289 


0.020 


99.222 
95.533 


B2 


Hi 
H2 


0.117 
0.183 


0.958 
0.697 


0.059 
0.348 


0.003 
0.042 


98.823 
98.784 


B3 


Hi 
H2 


0.229 
0.490 


0.593 
0.319 


0.779 
23.769 


0.103 
0.498 


98.408 
75.339 


B4 


Hi 
H2 


0.387 
1.066 


0.393 
0.166 


12.217 
36.597 


0.396 
10.313 


86.904 
52.879 


B5 


Hi 
H2 


0.921 
1.066 


0.188 
0.166 


36.500 
36.597 


6.787 
10.31 


56.475 
52.879 



Table 2: Branching ratios (in percent) and total widths (in units of GeV) for the 
Higgs bosons, H.i{i = 1,2), for the benchmark points of Tabled! Branching ratios 
that are not displayed, account for less than 0.4%. 

4.1 Zi7-production 

The first search channel for an invisibly decaying Higgs boson at the LHC considered here 
is the associated production of a Z and a Higgs boson, where the Z decays leptonically. 
This ensures that a corresponding event can be triggered. The backgrounds to this process 
include ZZ, WW, WZ and Z production with corresponding decays, and fully leptonic ti 
productioiij. It should be noted here that in principle some information on the rates can be 
obtained directly from data: for ZZ pairs, final states with four leptons may be reweighted 
with the corresponding Z — > z/z/ branching ratio, in the WW case, different sign, different 
lepton pairs may be invoked. For the WZ background, it may be possible to extrapolate 
from events where three leptons are seen to those where one lepton is lost, i.e. either outside 
the detector acceptance or undetected. For top-pair production, semi-leptonic events may 
help. 

All processes have been simulated with Sherpa [44] in the following setup: In order to 
correctly model hard parton radiation Sherpa employs the multijet matrix element-parton 
shower merging procedure of [45]. Therefore, for all processes discussed here and in the 
next section, matrix elements with at least one and in most cases two additional jets have 
been added to the simulation. This ensures that the simulation correctly describes the 
important high-p^ tails of various distributions. However, all cross sections quoted are, in 



Note that, in all processes, off-shell effects, Z-7 interference etc. are fully included in the simulation. 
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principle, obtained at leading-order accuracy, with no i^-factors added to them. CTEQ6L 
parton distubution functions are used with as{Mz) = 0.118 [46]. as is computed at two- 
loop accuracy. All scales are set according to the merging prescription of [45]. Jets have 
been defined in all cases through the krp algorithm [47]. The CKM matrix has been choosen 
to be diagonal. 

We have simulated and analysed events with electrons in the final state; mostly identical 
numbers would have been obtained if we had specialised for muon pairs instead. Obviously, 
this difference would be of great importance if detector effects had been included as welLJ. 
However it should suffice to state that we quote final results for leptons i = e, fi. We also 
omitted all effects due to the underlying event because of the large uncertainties related to 
its modelling and the rather small impact it has on the observables we discuss. 

The selection cuts listed in Ref. [28] have been applied. Thus we require: 

1. one lepton pair of the same kind with opposite charges, where each lepton individually 
satisfies pT,e > 15 GeV and |?7^| < 2.5; 

2. \Mu-Mz\ < 10 GeV; 

3. ^T > 100 GeV; 

4. a veto on jets with pt > 20 GeV, \ri\ < 4.9; 

5. a veto on b-jets with px > 15 GeV, \ri\ < 4.9; 



6. m^ > 200 GeV, where mx = y 2j9y|^r(l — cos0). 
Additionally, we impose: 

6. ARu < 1.75; 

7. priiif^T) < 60 GeV. 

For the various backgrounds listed above, cross sections before and after these additional 
selection cuts are listed in Table [31 Generation cross sections, selection cut efficiencies and 
the resulting selection cross sections for the signal in the different benchmark scenarios 
are given in Table HI It should be stressed again that all cross sections quoted have been 
obtained at leading order accuracy. 



^ We refrained from including full detector simulations, or any Gaussian smearing or electron-jet conver- 
sion "by hand" and concentrated on an analysis at the hadron level, including all effects of fragmentation, 
hadron decays, final state QED bremsstrahlung etc.. 
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zz 


W^Z 


w+w- 


ti 


Z 


^foT m 


164 


1.17-103 


1.01-10^ 


7.44-10^ 


1.81-10^ 


i+i- only 

\mu -Mz < 10 GeV 

^T > 100 GeV 

jet veto 

niT > 200 GeV 

ARu < 1.75, pt{U,:^t) < 60 GeV 


2.00-10-1 
1.87-10-1 
3.69-10-2 
1.64-10-2 
1.54-10-2 
1.23-10-2 


1.10-10-1 
9.17-10-2 
1.10-10-2 
2.13-10-3 
1.95-10-3 
1.50-10-3 


6.59-10-2 
8.92-10-3 

5.91-10-^ 
3.53-10-5 
2.74-10-5 
2.23-10-9 


8.40-10-2 
1.09-10-2 
2.41-10-3 
2.00-10-6 
1.19-10-9 
1.55-10-1° 


1.41-10-1 
1.25-10-1 
1.94-10-^ 


O-eflf [fb] 


2.02 


1.75 


2.25-10-5 


1.15-10-5 


- 



Table 3: Generation characteristics for the background processes to the ZH- 
channel. In all cases we included all leptonic decay modes: In the ZZ case, 
therefore the final state included a lepton and a neutrino pair, in the WZ case, 
we included a lepton pair from the Z and a lepton-neutrino pair from the W, the 
WW channel was supposed to decay fully leptonically in all possible combinations, 
for the top pairs we assumed purely leptonic decays, and for the Z a leptonic final 
state (no neutrinos) was demanded. 





Bi 


B2 


Bs 


B, 


55 


o-tot [fb] 


280 


114.6 


53.0 


29.0 


13.6 


i+i' only 

\mu -Mz < 10 GeV 

^T > 100 GeV 

jet veto 

niT > 200 GeV 

ARu < 1.75, pt{U,:^t) < 60 GeV 


1.75-10-1 
1.62-10-1 
3.12-10-2 
3.00-10-2 
2.88-10-2 
2.55-10-2 


1.98-10-1 
1.84-10-1 
6.07-10-2 
5.66-10-2 
5.49-10-2 
4.93-10-2 


2.25-10-1 
2.10-10-1 
8.91-10-2 
7.85-10-2 
7.64-10-2 
6.94-10-2 


2.40-10-1 
2.23-10-1 
1.08-10-1 
9.29-10-2 
9.08-10-2 
8.35-10-2 


2.24-10-1 
2.08-10-1 
1.12-10-1 
1.08-10-1 
1.06-10-1 
9.85-10-2 


O-eflf [fb] 


7.15 


5.65 


3.68 


2.42 


1.34 



Table 4: Generation characteristics for the signal processes in the ZH-ch.a.nnel. 
In each case we assumed all leptonic decay channels for the Z boson. 
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The numbers from both Tables [3] and H] suggest that the two most dangerous back- 
grounds to the ZH signal are ZZ and WZ production, with corresponding decays. Fol- 
lowing our discussion above, however, it seems that the total cross sections and distribu- 
tions related to these backgrounds can be directly extracted from data in the ZZ case or 
probably well extrapolated from measurements. After cuts we find that the backgrounds 
together account for roughly 8 fb, leaving us with signal-to-background ratios of the order 
of S/ B K, 1/8 up to 1. We therefore conclude that it should be possible to find the signal 
in all five benchmark scenarios. However, we would like to stress here that more conclusive 
numbers can be obtained after a simulation at detector level only. 

Such detector-level studies for an invisibly decaying Higgs boson have been for the 
ATLAS experiment [28] found signal-to-background ratios reaching up to 1/4. Although 
this is of the same order of magnitude as our results, there are several differences: First of all, 
in our simulation the Sherpa Monte Carlo event generator with multijet merging was used 
for both signal and background events, while the ATLAS study employed the Pythia [48] 
event generator for the backgrounds and the program h2hv [49] for the signal. While 
Sherpa and Pythia are formally of the same accuracy there are a number of differences, 
like Sherpa multijet merging leading to an improved treatment of hard QCD radiation, 
and the full inclusion of spin correlations in SHERPA, which are not present in PYTHIA. 
This may have lead to a better separation of signal and background in SHERPA. On the 
other hand, in ATLAS' simulation the HVV couplings where assumed to have exactly the 
same strength as in the SM - which is not true for our analysis, where these couplings 
are reduced due to mixing effects. In addition, a 100% branching ratio of Higgs boson to 
invisible was assumed for the ATLAS simulation, again in contrast with our simulation, 
where the relevant branching ratio ranged between roughly 50% up to 100%. These two 
facets of the study, of course, enhance the signal-to-background ratio in the ATLAS study. 
Of course, there are further differences, like the missing underlying event in SHERPA, which 
has been included in the ATLAS study, like slightly different selection cuts, like a different 
choice of PDF (CTEQ5L in ATLAS, CTEQ6L in our study) and, most importantly, like 
the inclusion of detector effects through their fast detector simulation ATLFAST [50] in 
the ATLAS study that are totally absent in our case. To summarize: However different 
in detail the studies are, it is reassuring to see that in all cases this seems to be a feasible 
channel, at least at accumulated higher luminosities. 

In addition to the findings above, cf. Tables [3] and H] we have identified two further 
distributions that may be worthwhile to study in the ZH channel: 

• The total transverse momentum of the leptons and /^r, i-e. the total transverse mo- 

20 
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Figure 6: The Pt{g^g~ iPt) distribution for the signal in benchmark scenario Bl 
and the individual backgrounds. The left panel displays individual distributions 
while the right panel show the sum of backgrounds and signal, starting from the 
lowest significant background. 



mentum of the H and Z candidates (see Fig.[6j). This observable shows a significantly 
different behaviour between the signals and the backgrounds, where the signal tends 
to be more strongly peaked towards small values. 

• The azimuthal angle between 'pT and the momentum of the lepton pair (see Fig. 
[7]). Here the signal tends towards a more back-to-back configuration of the Z and H 
candidate. Seemingly, there is a significantly higher QCD activity in the backgrounds 
than in the signal, providing more jets for the Zif -candidate pair to recoil against in 
the backgrounds. 

These findings may help to further improve the signal-to-background ratio. 

4.2 Vector Boson Fusion 

The other production channel we consider for invisibly decaying Higgs bosons at the LHC 
is vector boson fusion (VBF). As the name suggests, in this process the Higgs boson is 
produced through the fusion of two vector bosons emitted by quarks inside the protons, 
which typically carry comparably large momentum fractions of the protons. Therefore, 
at leading order (tree-level) there is no colour exchange between the two protons, and it 
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Figure 7: The A0(e+e~,^T) distribution for the signal in benchmark scenario Bl 
and backgrounds. The left panel displays individual distributions while the right 
panel show the sum of backgrounds and signal, starting from the lowest significant 
background. 

can be expected that the central rapidity region remains to a large extent empty apart 
from the decay products of the produced system. The quarks on the other hand will 
be deflected, typically by transverse momenta of roughly half the mass of the produced 
system. This gives rise to two hard jets, which, due to the invisible nature of the Higgs 
boson, are essentially the triggers in this analysis. The main background processes to be 
taken into account are the production oi Z or W bosons in association with two jets, 
which can originate either from QCD or through electroweak interactions, thus mimicking 
the topology of the VBF signal. In addition, top-pair production with subsequent semi- 
leptonic decays must be considered. Similar to the case of W production, the lepton is then 
lost. Again, it is worth noting that it should be possible to extract information concerning 
the total rates of these backgrounds, even after selection cuts, directly from data. This is 
possible either by reweighting leptonic Z decays to those into neutrino pairs, or, with a 
somewhat larger error, by extrapolating the modes where the individual lepton is seen (in 
ly+jets or semileptonic top-pairs) into those regions where the lepton is lost. This is in 
analogy to the case discussed above. We employ the basic cuts listed in Ref. [30], i.e. we 
require: 
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Figure 8: Tagjet invariant mass distribution for benchmark point Bl and the Z 
and W backgrounds. 

1. two tagging jets with 

(a) ptj > 40 GeV, |?7j| < 5, 

(b) \r]j^ -r/jj >5, r/j, ■ r]j^ < 0, 

(c) mj^j2 > 1700 GeV , 

(d) A0j^j2 = |0j^ -0j2| < 1, 

2. missing transverse momentum, 'f)T > 100 GeV; 

3. no identified lepton, i.e. no lepton with 
p"/ >5 ,6 GeV in \r]i\ < 2.5, 

4. a central jet veto, i.e. no jets with 

Pt > 20 GeV, m.m{r]j-^,r]j^} < i] < max{?7j^,77j2}. 

Additionally we impose: 



5- l^sl = \Vh - 1 iVh + Vj2)\ > 1-5, 
6. A0,,,,3, A0,„,3 < 1.25. 
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Figure 9: Left panel: Ht distribution for benchmark point Bl and Z background. 
Right panel: r^g- distribution for benchmark point Bl and Z background. 

The choice of the cut on the invariant tagging jet mass of mjj > 1700 GeV is motivated 
by the corresponding invariant mass spectrum shown in Fig. [HI We observe that the signal 
distribution crosses the background at rrijj ^ 1700 GeV. Of course this statement sensitively 
depends on the model parameters chosen; however, the common feature of all scenarios is 
that the higher the invariant mass cut, the better the signal-to-bachground ratio. This is 
due to the fact that in a large fraction of background events the two tagging jets originate 
from QCD or the decay of weak gauge bosons. 

After the above cuts the possibilities to check for the signal topology are limited in the 
VBF channel. Possible objects to be identified experimentally are the tagging jets, "^t and 
an eventually arising soft third jet. Therefore most observables show the same behaviour 
for signal and background, which is exemplified in the left panel of Fig. [9l showing the 
i^T-distribution for the signal at benchmark point Bl and the Z background. In the right 
panel of Fig. [9] we show for the same scenario the t]^ distribution. It is clearly seen that for 
the background the third jet tends to be more central between the tagging jets, while for 
the signal it is rather forward or backward. This motivates the first of the additional cuts 
above. 

For the various backgrounds listed above, cross sections before and after additional 
selection cuts, and the number of generated events are listed in Table [5l Signal cross 
sections before and after additional selection cuts are listed in Table [61 Putting together 
numbers, we again find appreciable signal-to-background ratios between more than 1/3 
up to nearly 1 for all the benchmark points in the model. However, this finding has to be 
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Z+jets 
(QCD+EW) 


PF+jets 
(QCD+EW) 


tt 


^foT [nb] 


9.41 


51.8 


0.145 


tagging jets 
rrijj >1700 GeV 
;^T > 100 GeV 
lepton veto 

central jet veto 
r]*\ > 1.5 


1.80 
3.49 
2.64 
2.63 
4.03 
1.54 
1.37 
1.14 


10-4 
10-5 
10-5 
10-5 
10-6 

10-6 
10-6 
10-6 


7.44 
1.64 
9.73 
2.84 
9.87 
2.18 
1.95 
1.44 


10-5 

10-5 
10-6 
10-6 
10-7 
10-7 
10-7 
10-7 


1.62 
4.44 
3.32 
1.28 
2.79 
1.44 
6.70 
4.29 


10-3 
10-4 
10-4 

10-5 
10-5 
10-6 
10-7 
10-7 


o-cff [fb] 


10.7 


7.45 


0.0621 



Table 5: Generation characteristics for the background processes to the VBF- 
channel. Here, the Z boson decays to neutrinos, whereas the W boson decays to 
any lepton-neutrino pair. For the top-pairs, semileptonic decays only have been 
considered. 





5i 


B2 


Bs 


B, 


B, 


f^toT [pb] 


5.46 


4.46 


2.99 


2.06 


1.32 


tagging jets 


4.38 


10-2 


5.59 


10-2 


6.80 


10-1 


7.54 


10-2 


8.12 


10-1 


mjj >1700 GeV 


1.69 


10-2 


2.20 


10-2 


2.74 


10-2 


3.07 


10-2 


3.36 


10-2 


I^T > 100 GeV 


1.46 


10-2 


1.90 


10-2 


2.37 


10-2 


2.67 


10-2 


2.92 


10-2 


lepton veto 


1.46 


10-2 


1.90 


10-2 


2.37 


10-2 


2.66 


10-2 


2.92 


10-2 


^^nn <1 


5.76 


10-3 


7.65 


10-3 


9.46 


10-3 


1.06 


10-2 


1.18 


10-2 


central jet veto 


3.42 


10-3 


4.33 


10-3 


5.35 


10-3 


6.06 


10-3 


6.64 


10-3 


r]*\ > 1.5 


3.40 


10-3 


4.31 


10-3 


5.32 


10-3 


6.03 


10-3 


6.60 


10-3 


A0j,j3, A0J2J3 


3.11 


10-3 


3.92 


10-3 


4.81 


10-3 


5.45 


10-3 


5.97 


10-3 


o-cflf [fb] 


17.0 


17.5 


14.4 


11.2 


7.9 



Table 6: Generation characteristics for the signal processes in the VBF-channel, 
for the different benchmark scenarios. 
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taken with more than a pinch of sah: first of all, similar to the ZH channel, we included all 
effects due to fragmentation, hadron decays, QED bremsstrahlung etc., and we typically 
added at least one further jet for a better modelling of additional hard QCD radiation. 
We did not, however, include the effects of the underlying event, which here could play 
a significant role in filling the rapidity gap between the two taging jets, and thus lead to 
a corresponding reduction in the effective cross section after selection cuts. In addition 
we did not include diagrams where the Higgs boson is produced through an effective ggH 
coupling, mediated by heavy quarks. Although in principle the cross section for this mode 
is large, we note that previous work in the framework of the Standard Model suggests that 
the typical VBF cuts render this contribution insignificant [51-54]. Also, again, we did not 
simulate events at the detector-level which could further modify our findings. 

However, again our results are in qualitative agreement with results of such a simulation 
at the detector level, which has been performed for the ATLAS experiment [30]. The results 
of this study were obtained using a fast detector simulation, and they are quite encouraging, 
too. Although in quantitative agreement, there are several differences: Again, the first one 
lies in the choice of the evenet generator. ATLAS chose Pythia to compute both signal and 
SM backgrounds at leading order, while we employed Sherpa. In the ATLAS simulation SM 
coupling strength for the HVV couplings has been assumed with a 100% branching fraction 
of the the Higgs decay to invisible, while in our study the HVV coupling is shielded through 
the mixing of the scalars, and the relevant branching ratio ranges between 0.5 and 1. While 
in ATLAS' Pythia simulation the effect of hard QCD radiation is typically accounted for by 
the parton shower, Sherpa uses exact matrix element, leading to a significantly increased 
jet activity. Also, SHERPA naturally includes spin correlations, and VBF-like background 
topologies are also taken care off, which have been missed in the ATLAS simulation. These 
effects, together, would typically reduce the signal-to-background ratio in our simulation 
with respect to the ATLAS study. On the other hand, the effect of the underlying event as 
well as the fast detector simulation, both included in the ATLAS analysis but ignored by 
us, may have the opposite implications on the visibility of the signal. Finally, it is worth 
stressing that we have also chosen different optimization cuts, in particular cuts 5 and 6, 
to enhance the signal over the background. 

Nevertheless, to summarize, we again find that the prospects of finding an invisibly 
decaying Higgs boson at the LHC are much better than naively anticipated, and the two 
channels considered here may very well play a significant role in the phenomenology of 
non-standard scalar sectors. 
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5 Non-Abelian Phantom Sector 

So far only a Gp = U{1) group theoretic phantom sector has been considered. The obvious 
question to be asked is how the Higgs boson observabihty will be affected in the case of 
non-Abelian extensions of the phantom sector (like Gp = SU{N)). This will briefly be 
discussed in this section. As an overall result, in general, such extensions typically result 
in further suppression of the Higgs boson visible event rates, 7?.^. Furthermore, in the 
case of more involved representations or multiple vector representations of Gp the "Higgs 
—>■ invisible" signal is decreased to a non-detectable rate. Some examples supporting this 
result will be presented in the following. 

— * 

Consider for instance a Gp = SU{N) vector representation of scalar phantom fields, $. 
Then SU{N) is spontaneously broken down to SU{N — 1) with 2N — 1 physical NGBs and 
one physical SM-singlet scalar field that eventually mixes with the SU{2)l Higgs field. It 
is a textbook exercise to prove that eq. (11. 5p in such a framework becomes 

2 
777 

Ant = --^Oi2H,{x)J''ix)J''ix) with a = l...{2N-l). (5.8) 

/ a 

This suggests that the Higgs boson decay width broadens compared to the Gp = f/(l) 
case. The visible Higgs boson event rates (there are still two physical states) read 

-1 



nl 



nl 



9 , / 2N - 1 m?r, „ „ 

1 + tan^ ^ 1 + ^ tan^ 6 tan^ 8 

V 12 mi 



(1 + cot^ 9) (l + =^—^ ^ cot^ e tan^ /3 
\ 12 mi 



Hence increasing the rank of the phantom gauge group results in a (1/A^ for large A^) 
decrease in visible Higgs boson rates. Searching for "Higgs -^ invisible" is therefore vital. 
Note also that increasing the rank of the phantom symmetry group does not necessarily 
imply different "Higgs -^ invisible" rates. In fact, in the above example we still have two 
physical scalars in the spectrum for which the equation 7i + 72 ~ 1 is valid, similarly to 
the Gp = f/(l) case. 

It may also be the case that additional physical Higgs bosons fragment the "Higgs — > 
invisible" rate into many small pieces such that any detection at the LHC seems completely 
impossible. This case can be illustrated with the following example: consider Gp = SU{3) 
broken by 2 sets of vector representations down to the null group. We start with 12 degrees 
of freedom, out of which 8 become NGBs and the other 4 become massive scalar fields. 
These 4 fields will mix with the one SU{2)l Higgs field through the (5 x 5) matrix O 
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forming 5 physical Higgs-boson eigenstates. In this case, due to the unitarity of the matrix 
O we have X]j=i ^^ ~ 1; which allows for T^ < 0.25. Such a "Higgs — >• invisible" rate 
is most probably beyond reach of discovery (or exclusion) at the LHC [28, 30] - a truly 
nightmarish scenario! 

6 Additional Remarks 

It should be emphasized that in the scenario considered in this article, invisible Higgs 
boson phenomenology, small neutrino masses and the correct baryon asymmetry (see also 
ref. [13]) are all obtained without fine-tuning coupling constants. All scalars have masses 
at the EW scale (tan/3 ~ 1) and so there are no ultra-heavy scalars to destabilize this 
hierarchy. However, the model does not include gravity nor does it contain a mechanism 
or theoretical explanation as to why a <^ Mpianck- Although the SM hierarchy problem 
is not solved in this model the question here is somewhat different: Why is the phantom 
sector symmetry broken at the EW scale ? We cannot provide a non-common {i.e., non- 
super symmetric) answer to this question, and refer to [55,56]. 

Instead of a theory with one global symmetry, one could imagine a theory where sev- 
eral symmetries were gauged (or left un-gauged), absorbing the NGBs into massive gauge 
bosons through the Higgs mechanism. This is an absolutely viable option, although the re- 
quirement of anomaly cancellation would result in model dependencies. Such models have 
been proposed before and studied in some detail in the recent literature [57]. Generally 
speaking, these models lead to phenomenology that includes the (observable) decays of the 
extra gauge bosons, with all constraints on their masses etc.. 

Recently there has been renewed interest in the possibilities offered by extending the 
Standard Model with a real scalar singlet [58]. Depending on the symmetries of the model 
it is possible to provide a candidate for the cold dark matter in the universe (extra discrete 
symmetries needed) [59] , and it is possible to provide a strong first-order electroweak 
phase transition suitable for electroweak baryogenesis [19]. It should be noted in the latter 
case that an additional source of CP-violation would be necessary to provide a complete 
mechanism for baryogenesis. 

Models with broken discrete symmetries provide another possible way of avoiding in- 
visible decays of Higgs boson(s). Clearly, spontaneous breaking of such symmetries does 
not lead to NGBs, making the Higgs boson signatures more visible. There are, however, so 
many possibilities of such groups that a particular choice renders this idea less appealing and 
convincing. Spontaneously broken discrete symmetries may also, in some cases, produce 
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unwanted cosmological relics such as domain-walls, potentially placing severe constraints 
on this class of model. 

7 Conclusions 

Physical NGBs arise when continuous global symmetries are spontaneously broken. Such 
broken symmetries may be related to the smallness of neutrino masses or the patterns of 
mixing angles (in the case of familons). In this article we show that the role of NGBs 
in Higgs boson phenomenology is very important; they lead to the dilution and potential 
invisibility of the expected SM signal. Working with approximate analytic formulae we 
first identified regions of parameter space [eq. (11. 6p ] where Higgs boson phenomenology 
is challenging both for past LEP data and for the future LHC experiments, and secondly 
implemented the model in Sherpa, ready for further analysis by experimenters when real 
LHC data arrive. 

Our study shows that LEP excludes the minimal phantom sector case where both 
Higgs bosons have masses mu < 85 GeV irrespective of their decay modes. However, 
experimentally allowed scenarios exist where one Higgs boson mass is much lower than 
the SM Higgs boson exclusion limit, rriHi = 68 GeV, and the other is just at this limit, 
771^2 = 114 GeV. 

In light of the nightmarish potential of this scenario, Monte-Carlo simulation studies of 
invisible Higgs boson searches at the LHC are performed. Two search channels are looked 
at in detail; the associated production of a Z and a Higgs boson, and the production of a 
Higgs boson in weak vector boson fusion. For ZH associated production, it is found that 
in each of 5 benchmark scenarios, the invisible Higgs boson should be found at the LHC, 
with signal-to-background ratios of order S/B ^ 1/8 to 1. Scope for improving this ratio 
is also found by looking either at the distribution of the total transverse momentum of the 
leptons and the ^t, or at the distribution of the azimuthal angle between the '^t and the 
momentum of the lepton pair. Fairly good signal-to-background ratios are also found in 
the vector boson fusion search channel. However in this case the effects of the underlying 
event, which was not included in simulations, may reduce the amount of signal passing the 
selection cuts. 

Although our MC analysis focuses on the case with an Abelian phantom sector sym- 
metry, we also examined cases with non- Abelian symmetries in the phantom sector using 
the analytic formulae provided in section [31 For the case Gp = SU{N) we found that the 
visibility of the Higgs bosons is reduced when we increase the rank of the SU{N) group 
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making the LHC searches to invisible a necessity. In addition, by choosing appropriate 
representations of the group for breaking the symmetry we may further dilute the Higgs 
boson to invisible signature, leading to a very difficult scenario indeed for the LHC. 

Regarding the hierarchy problem, the model at hand is not better or worse than the 
Standard Model. Any difference could be interpreted as shifting the problem to the phan- 
tom sector which sets the scale of the symmetry breaking. 
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Appendix A ^(1) Phantom Model Feynman Rules 

In this appendix we present Feynman rules for the Higgs sector of Gsm x {Gp = f/(l)p} 
that are relevant for Higgs phenomenology at LEP and the LHC. Feynman rules for the 
trilinear couplings HiJJ, HiH2Hi, HiHiHi, W+W-Ri, ZZH^, and ffHi for i = 1,2 are 
shown in Fig. [101 For completeness in Fig. [HI also Feynman rules for the quadrilinear 
couplings HiHjJJ, HiHjHkHi, HiHjZZ, and HiHjW+W' are listed. 

M\Y and Mz are the masses of the W boson and Z boson, respectively and rrif is the 
fermion mass which can be either a quark or a lepton. The SU{2)l coupling constant is g2 
and 6w is the Weinberg mixing angle, v is the vacuum expectation value for the standard 
model SU{2)l Higgs doublet H. g^^ is the Minkowski spacetime metric (1, —1, —1, —1). 

The orthogonal mixing matrix, O, is 

^^/On 0.\ / COS. sm.y 

Here tan/3 = v/a with a = ($). niH^ and ttih^ denote the masses of the two Higgs bosons. 
Hi and H2, respectively. 
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Figure 10: Trilinear couplings. 
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Figure 11: Quadrilinear couplings. 
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